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Abstract. We derive the conditional probability distri¬ 
bution /pk(t> | v) of the peculiar velocity around a peak 
of given overdensity v in a Gaussian density field. This 
distribution characterizes the shear field around local, iso¬ 
lated density peaks. We use the moments of the velocity 
distribution to study the spatial distribution of the escape 
(jrsjlocity from protostructures forming in a standard, flat, 
©Nile-invariant CDM model. In the case of high (v ~ 2) 
i-psiaks with filtering scales Rf ~ 2/) _1 Mpc, we find that 
distances between 2 4- 4 Mpc h 1 (or equivalently k, 
1,4 -2 Rf) a sensitive fraction of the matter enclosed within 
a! i extrapolated (up to the present epoch) volume contain¬ 
ing a total gravitating mass M ft* 2.725 x 10 15 /) _1 Mq (of 
JjBe same order as that derived from gravitational lensing 
Reservations, see Tyson and Fisher 1995) is gravitation- 
puly unbound to the peak. The fraction of unbound matter 
(Tctries between approximately between 0.05 — 0.2 depend- 
Ufi^ on the central height v of the density peak, on the 
^hering radius and on the adopted normalization for the 
flu ctuation spectrum. We show how these results support 
Gkl previous conclusions regarding the delay of the col- 
}y se induced by the presence of small-scale substructure 
fej clusters of galaxies formed in standard CDM cosmolog- 
S03d models. 
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1. Introduction 

In most of the contemporary, widely studied cosmologi¬ 
cal scenarios the formation of cosmic structures is traced 
back to the evolution of primordial density fluctuations. In 
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particular, within the hierarchical clustering scenarios like 
those envisaged by Cold Dark Matter (hereafter CDM) 
models, local peaks of the density field are regarded as 
sites which could have engendered present-day observed 
structures like galaxies and clusters of galaxies (Kaiser 
1986; Bardeen et al. 1986, hereafter BBIvS). The struc¬ 
ture of the density peaks arising from primordial Gaus¬ 
sian perturbations has been studied by many authors (e.g. 
Peacock A Heavens 1985; BBIvS; Gonzalez A Martinez 
1995), with particular attention to their mass distribution. 
The structure of the velocity field around a peak, however, 
has received comparatively lesser attention, although its 
knowledge can be important for the interpretation of the 
observed bulk flows (Davis A Peebles 1983; Gorski 1988; 
Groth et al. 1989) and for a quantitative understanding of 
the collapse and virialization of protostructures (Weygaert. 
A Babul 1994). The first, point is particularly important 
when one compares the predictions of cosmological mod¬ 
els with peculiar velocity data. Recently Mo et: al. (1993) 
and Marzke et al. (1995) have shown that the discrepancy 
beteween the standard CDM predicted galaxy pairwyse 
velocity dispersion and the observed value is greatly re¬ 
duced when the contribution from galaxies in rich Abell 
clusters is subtracted from the samples. It remains to be 
proved that the predicted contribution from these clusters 
in standard CDM models is of the right order of magni¬ 
tude to explain the discrepancy. In this paper we will show 
that this is indeed true, although with a rather large vari¬ 
ance. 

A knowledge of the velocity distribution function around 
a density peak can be of importance also to give quantita¬ 
tive information on the amount od substructure around a 
proto-object. Indeed, in some recent papers (Antonuccio- 
Delogu, 1992; Antonuccio-Delogu A Colafrancesco 1994; 
Colafrancesco et al. 1995) we have shown that dy¬ 
namical friction induced by the presence of small-scale- 
substructure (hereafter sss, a natural byproduct of any 


delay of the collapse of the outermost shells of a proto¬ 
cluster. However, we also remarked in those papers that 
only that fraction of the sss made of subpeaks not gravi¬ 
tationally bound to the perturbation should be considered 
as part of the “background” and only its contributions 
should enter into the calculation of the dynamical friction 
coefficient. The conservative estimate we have adopted for 
this fraction was ~ 0.1. In fact, this fraction can be easily 
computed once the velocity distribution function is known. 
As we will see, the calculations we report in this paper will 
prove that our estimate was indeed a conservative one. 
The plan of the paper is the following. In Sec. 2 we will 
derive the probability distribution of the velocity around 
peaks of given central overdensity v, using standard meth¬ 
ods from Gaussian random Helds’ theory (as sketched, e.g. 
in the Appendix D of BBKS). This distribution is gener¬ 
ally anisotropic, and its knowledge fully specifies the ve¬ 
locity field around the peak. In Sec. 3 we will then use 
this knowledge to determine the fraction of gravitation¬ 
ally unbound matter. Finally, in Sec. 4 we will present the 
conclusions and some hints for future work. 


2. Velocity probability distribution around a den¬ 
sity peak 

2.1. Conditional probability distributions 

Our target is to calculate the peculiar velocity probabil¬ 
ity distribution f p k(v | v) around a peak of given central 
overdensity v = bo/cr. The peculiar velocity v is a random 
variable given (in the linear regime) by: 


The above probability distribution is a conditional mul¬ 
tivariate probability (under the hypothesys that all the 
above variables are Gaussian distributed), and can be 
computed using standard theorems. Let us define Ya = 
(v,r)), Yb = (v,x). If Y”a,Y”b are Gaussian distributed, 
then the multivariate probability distribution for the vec¬ 
tor Yb given that the vector Ya has a given value, is a 
Gaussian with mean: 

(Y b |V A ) = (VB®VA)(VA®V A )- 1 vi (2) 

and covariance matrix: 


(AY b ® AY B | Y a ) = 

(Y b ® Y b ) - (Y B ® Y A )(Y A ® Ya)“ 1 (Y a ® Y B ) (3) 


(see BBKS, Appendix D). Here AY B = Y B - (Y B | Y A ), 
and the symbol (g) denotes a matrix product. 

The explicit calculation of these matrices is not difficult 
but a little tedious and can be easily done once some av¬ 
erages like (vv) are computed (to make easier the job to 
the interested reader we calculate all the relevant non-zero 
averages in the Appendix A). It turns out that the vector 
AYb is given by: 


AY B( i) 


( vi - v(v iv ) - JrEfci mivirii) ^ 
v 2 - v{v 2 v) - Jr EL 
V 3 - v{v 3 v) - 4r Ef=l Vi^i) 


( 4 ) 
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x — v{vx 
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and the final result for maxima (rj = 0, x > 0) reads: 


A standard theorem from the theory of gaussian random 
fields ensures that, being linearly dependent on a gaus¬ 
sian random variable (the overdensity 6(r)), the peculiar 
velocity will also be randomly gaussian distributed. 

We consider here relatively isolated, high v density peaks 
which are nearly spherically symmetric. The analogous 
probability distribution around two peaks of the den¬ 
sity field has been recently considered by Regos and Sza- 
lay (1995). 

We will first compute a more general distribution: f(v, x \ 
v, rf), i.e. the conditional probability for velocity and cen¬ 
tral curvature (x = — V 2 (5/<72 |o) f° r peaks of given v 
whose central shape is described by the gradient of the 
field, r) =< V6(r) >, where 6(r) is the density field found 
around the peak (see e.g. BBKS). 

We will then average this probability distribution over 
maxima (x > 0) imposing the condition r) = 0, char¬ 
acterizing points of the random process that are extrema. 
Note that f(v,x \ v,rf) denotes the probability distribu¬ 
tion of the 4-vector (v, x) given that the quantities (v, rf) 


f(v,x \i/,n = 0) 


1 

( 2 *) 2 (nLiM 1/2 


• exp 


E 





(5) 


In this latter equation we have introduced the eigenvalues 
A i,i = 1, . .4 of the covariance matrix M = (AY b <X> AY a) 
whose elements are: 


Mi, 


= n(r) - ! { Vi v } 2 + -j J2( ViT lj) E * = !> ••• 3 


j =i 


M 44 — 


tv 2 ) - (vx) 2 


Mi. 


= n(r) - E(r) - l (1>iv){ Vj v) + -3 
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Mi4 = 


i,j = l,...3,i^i 

(vix) - (viv)(vx), * = 1 ,...3 


( 6 ) 



functions II(r), £(?•) are the parallel and transverse com¬ 
ponents of the velocity tensor for the background field: 

(v(r')v(r' + r)) = [II(r) - £(?•)] r a r l3 + E(r)6 al3 (7) 


These functions enter our formula because of the rela¬ 
tion (i>iVj)(r) = (vi(r')vj(r' + r )} holding in an homo¬ 
geneous background. 

Formulas for II(r) and £(?•) are given in Gorski (1988). 
From a computational point of view we find more con¬ 
venient. to adopt, the equivalent, formulae given by Pee¬ 
bles (1993), [ his Eqs. (21.72) and (21.74)]: 
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with J n (r) and K n (r) given by: 


J n (r) = / dr£(r)i 
Jo 


n — 1 
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Fig. 1. Eigenvalues of the velocity correlation tensor for a stan¬ 
dard CDM model with bg = 1.1, a filtering radius Rf = 5 Mpc 
and v = 1.2. 


From BBIvS [1986, see their eq. (A3)] we have: (vx) = j. 
The other averages appearing in the various terms of the 
matrix are given by: 

n ( dr'r' 2 £(r'), ( 10 ) 

r (To Jo 


{ v iVj) =-^— r-ij-j yX(r)~ J ar $(r j j , 


{'Vix) 


H 0 Cl° 0 6 d£ 

- r i-r- 

V(J2 dr 




( 12 ) 


(see the Appendix A for more details on the derivation of 
these quantities). 

From the distribution function given in Eq. 5 we can ob¬ 
tain i In average velocity distribution around peaks of ar¬ 
bitrary height, by imposing the condition = 0 and inte¬ 
grating over all values a? > 0. Eventually we obtain: 


/ P k(v | v) = 


[a/At 


(IlLi W 1 * 

where we have defined: 


exp 


"=—" ,c 




2 = 1 
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(13) 



Fig. 2. S ame as Fig. 1 for v = 2 and for Rf = ... such that the 
total mass at. a given distance of lOMpc is M PS 2.725 x 10 15 Mg. 

From Eq. (13) the physical meaning of the eigenvalues 
A i,i = 1, ...3 is clear: they are: the components of the pe¬ 
culiar velocity dispersions along the parallel (Ai) and per¬ 
pendicular (A 2 3 ) w.r.t. the radial direction of a. density 
peak. In Figs. 1 and 2 we plot, the profiles of these quanti¬ 
ties for two configurations having approximately the same 
total mass (M ps 2.725 x 10 1 5 A Mq ) at a. reference distance 
of 10 /) _1 Mpc (comoving units). This choice of mass and 




total cluster masses as derived from gravitational lensing 
observations in nearby clusters (Kneib ft al. 1995; Wax- 
man A Miralda-Escude 1995). For example, in Abell 1689 
(Tyson & Fischer 1995) a mass M ps 2.7 x 10 15 /) -1 Mq 
is found within Ra = 1.5 /) -1 Mpc. The collapse factor 
of a density shell reaching an overdensity 6f is given by 

x c = [l + 6f ] 1/,S (Ryden & Gunn 1987), and for 6f = 300 
we get x c ps 6 .6. Thus, the final radius of a shell at 
10/) -1 Mpc is approximately equal to one Abell radius, for 
values of the central overdensity typical of rich Abell clus¬ 
ters. High v, i.e. more centrally condensed density peaks 
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Fig. 3. S ame as Figure 1 for bs = 1.8. 

have a larger spatial gradient of the velocity dispersion, as 
one can intuitively expect. It is however interesting to note 
the considerable quantitative difference among the central 
values of the velocity dispersions: about. 100 Km/sec. This 
effect clearly suggests a broad variance in the magnitude 
of the expected bulk flows. 

2.2. Dispersion profiles 

The three components of the velocity dispersion attain 
nearly the same value in the central regions of the peaks, 
within the filtering radius (in Fig.l Rf = 5 /) -1 Mpc while 
in Fig.2 Rf ~ 2 h -1 Mpc). The velocity dispersion Ai in 
the radial direction decreases faster with the distance r 
from the center of the peak than the tangential velocity 
dispersions Ao and A 3 . In fact, at larger radii r, the ra¬ 
dial component is sensibly lower than the tangential com- 


of the clusters are predominantly non-radial. For objects 
with the same mass, the onset of non-radial motions hap¬ 
pens at shorter radii when the filtering radius is smaller 
and consequent.ely the central height of the peak is higher 
(see Fig.2). 

In 1 lie wry central regions of the peaks, r A 0.5 Rf, there 
is no net difference among the three components of the 
velocity dispersions and they are approximately constant 
meaning that this region is dominated by an isotropic ve¬ 
locity distribution. 

The region contained within the’ filtering radius does not 
show a net predominance of radial motions. As our consid¬ 
erations refer to protostructures still in the early phases 
of their non-linear evolution, this region can be probably 
considered as the one that will collapse in the final non¬ 
linear configuration through an infall process. 

Regions of the protost.ructures at r t> Rf contain mainly 
non-radial motions. In these outer regions the fate of 
the infalling material could be strongly dependent on the 
amount of tangential velocity dispersion relative to the 
radial, inward acceleration. This can be also seen by look¬ 
ing at the equation of motion of a spherically spherically 
symmetric mass distribution with density n(r): 

d . . d . o, 

— n (r)(v r ) = -n(r)— - —n(r)(v;) + 

„(,.) ( ”'>- 2(,a (14) 

r 

From the structure of the right hand side of equation (11), 
it is easy to see that high tangential velocity dispersions, 
{Vf) 2 ; 2 (t), 2 } may alter the infall patterns. 

Some recent 

numerical simulations (Evrard A Crone 1992), show that 
in the infall regime the velocity is primarily radial. These 
results were obtained using numerical codes with limited 
dynamical range; in fact, we believe that high resolution 
simulations are needed to determine the detailed struc¬ 
ture of the velocity Held in clusters of galaxies. But based 
on our results we are lead to conclude that in the early 
stages of the cluster evolution the velocity distribution is 
not prevalently radial, even tough the strong non-linear 
evolution tends to modify the velocity distribution of the 
infalling material. However, as clusters of galaxies - being 
the largest and youngest gravitationally bound systems in 
the universe - show a broad variance of dynamical states, 
the debate on the structure of their velocity distribution 
is still open. 

Finally, note that from Figs. 1-2 typical values for the 
velocity dispersions at 1 /) _1 Mpc are in the range 350 — 
450Km/sec, i.e. well within the 1 a values found by 
Marzke et al. (1995): 540 ± 180 km/sec. These quantities 
are not easily comparable to the pairwyse velocity corre¬ 
lations measured by the latter authors and by others (e.g. 
Davis A Peebles 1983), because pairs of galaxies can be 







Fig. 4. Spatial variation of F(\ v |> v vesc ) in the case of a 
top-hat perturbation of radius R = 1.5ft -1 Mpc. The three 
curves refer to values of the velocity dispersion, a v = 800 
km/s (dashed curve), <r v = 1000 km/s (continuous curve) and 
a v = 1500 km/s (dotted curve). Here ft = 0.5. 


found within region of different average overdensity. Even 
admitting that the galaxy pairs selected in a survey sam¬ 
ple fairly the (average) peculiar velocity Held of the regions 
to which they belong, one should expect anyhow a rather 
large intrinsic uncertainty at any scale because of the in¬ 
trinsic variance of the peculiar dispersion profiles. We will 
deal with these problems in a future paper (Antonuccio- 
Delogu et, al. 1995). 


3. Fraction of gravitationally unbound matter 


The previous considerations allow us to evaluate the frac¬ 
tion of matter that is gravitationally bound (in a statis¬ 
tical sense) to the protostructure. A simple criterion to 
define a gravitationally bound, infalling shell is to com¬ 
pare the velocity dispersion at distance r with the escape 
velocity v esc (r) = (2 | $(?•) I) 1 / 2 at the same distance 
from the center of the protostructure. 

A non zero fraction of unbound matter could be expected 
even in the case of a uniform system with a Maxwellian 
velocity distribution, f(v) oc exp(—v 2 /‘2a 2 ), characterized 
by a uniform velocity dispersion a 2 . In fact, in the case 


Mpc (the Abell radius), the fraction 


F( v > v esc ) 


I\v\>y^ d 3 v f^ 

f d 3 vf(v) 


( 15 ) 


takes values k, 0.08 -y 0.33 for a v = 800 -y 1500 km/s 
(see Fig. 4). This fraction rises up to values 0.4-P 0.7 at 
distances r ~ 5ft -1 Mpc for the same values of the velocity 
dispersion. 

For the case of a density peak, the fraction of unbound 
matter reads: 


F(v > v esc ) 


f d 3 vf pk (v,v) 


(16) 


We compute this quantity for a spherically symmetric den¬ 
sity peak in a standard CDM model. We have been able 
to reduce the integral in the numerator of Eq. 16 to a one 
dimensional integral which we evaluate numerically: the 
interested reader will find the details of this calculation in 
Appendix B. 

In Fig.5 we show the result in the case of peaks with dif¬ 
ferent heights subtending the same mass as given in the 
preceding section (M 10 ft-i.Mpe = 2.725 x 10 15 Mq). 



Fig. 5. Spatial variation of F(\ v |> c iiesc ). The thick solid 
curve is for bg = 1.6, the other two for bg = 1.1. The total 
mass for each curve is the same as in Figs. 1-3. 


For a relatively low peak v = 1.2 (which is quite fre¬ 
quent at cluster scales) we obtain that the fraction of un¬ 
bound clumps within Rf = 5ft -1 Mpc is F(v > v esc ) k, 








at ~ 2 Rf = 10/i -1 Mpc. This fraction lowers for higher 
v peaks, as the potential of these protostructure is deeper 
and hence the value of v esc is systematically higher. We 
found, in fact that, for a peak with v = 2 and Rf = 2/i -1 
Mpc the fraction F(v > v esc ) x 0.04-P 0.07 within Rj, 
and rises up to values 2 ; 0.1 at 2 Rf. 


4. Conclusions 

In this paper we have been mainly concerned with the 
structure of the velocity held that is found around the local 
maxima of the density fluctuation held at early times. We 
have derived analytically the velocity distribution func¬ 
tion around a density peak starting from the approach of 
BBKS, originally used to derive the density prohle of the 
peak in a Gaussian random held. 

From the analysis of the three components of the velocity 
correlation tensor, we found that peaks show the presence 
of a nearly isotropic region in the very central regions (at 
scales r x 0.5 Rf). Non-radial motions, instead, become 
dominant in the outer regions of the protostructures at 
r'ZRf. 

As the region of the protostructures that develops infall is 
found probably at R ’x, Rf , we are lead to the conclusion 
that in the early stages of the protostructure evolution 
there is not a coherent infall of matter during the collapse 
of the structure. This situation might involve a modifica¬ 
tion of the collapse dynamics of such structures, at least 
in their early phases. 

In the standard CDM scenario the flatness of the power 
spectrum at high wavenumbers produces an high fraction 
of clumps on small scales populating the external regions 
of the clusters. According to the previous results, these 
clumps can have high velocities (momenta), even larger 
than the average escape velocity from the structure. As 
the escape velocity is a measure of the cluster gravitational 
potential, the frequency of clumps that can be statistically 
of high momenta (those we call “unbound”) is larger in the 
outer regions of the protostructure. In fact, we evaluated 
that the fraction of unbound clumps within the filtering 
radius of the clusters is not negligible: we found values 
F(v > v esc ) ~ 0.05 -y 0.4 for quite common peaks, with 
central heights in the range v ~ 1.2 -P 2 . 

These findings indicate that there exist a large fraction 
of small mass clumps with statistically high momenta in 
the peripherical regions, r 2; Rf, of the protoclusters. 
During the collapse of the cluster, and also during the 
subsequent secondary infall, these high velocity clumps 
can exchange momentum with the collapsing or infalling 
clumps and generate a drag force (produced by dynami¬ 
cal friction) onto the collapsing material that acts to de¬ 
lay their collapse time scales (see Antonuccio-Delogu & 
Colafrancesco 1994, Colafrancesco et al. 1995). The set of 
high momentum clumps provides a Held of stochastic force 
generators that intervenes in modifying the equation of 


a delay of the collapse time particularly affecting the ex¬ 
ternal shells of the clusters. This reflects eventually in a 
collapse delay of the overall structure. Such considerations 
support and strengthen the conclusions of our previous 
work (Colafrancesco et al. 1995) concerning the dynam¬ 
ical effect of small-scale substructure on the collapse of 
clusters of galaxies. 
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A. Appendix: Velocity field correlations 

The ensemble-averaged vector products appearing in 
Eqs. ( 6 ) are linear combinations of the following averages: 

(x 2 ), {xv), (v 2 ), (vv), {vrj), (w) (Al) 

From BBKS, eq. (A3) we have that: 

{x 2 ) = 1, (xv) = 7 , (v 2 ) = 1 (A2) 

where: 7 = crl/a^cro and oy is the i-th momentum of the 
variance of the fluctuation field: 

Here P(k) is the filtered power spectrum. 

Remembering the definition: v = 8 (0) / (To and the defini¬ 
tion of correlation function: 


£(r) = («(r)«(0)) 


(A3) 


and adopting Eq. 1 for the peculiar velocity in the linear 
regime, we immediately find: 

drV 2 / 2 / (r') (A4) 

Let us now evaluate the quantity (vrj). The components 
of this tensor are given by: 

(viiy) = - g °^° n f drV ,2 (/(r')ry i ), 

r Jo 

where r)j are the components of the gradient of the density 
field evaluated at the center of the peak density profile: 


h 0 ng - 6 r- f 
r 3 (To Jo 


iy = 


d8(r) 


dr; 


(A5) 


Performing an integration by parts and recalling the defi¬ 
nition of correlation function (eq. A3) we obtain: 


, , H 0 n° 6 

{ V iVj} =- Z3 


H 0 Q°o 6 


dr 1 " 12 


r-WrRjJ 


dr'-' 2 3 


r = 0 

rj d£(r') 


H 0 n° 6 


"C(r) ~ 


.1 ^1 

dr'^(r') 


(A 6 ) 



Our purpose is to evaluate the integral in the numer¬ 
ator of Eq. 16 for the probability distribution given in 
Eq. 13. We make a change of coordinates in velocity space: 
(iq, i> 2 , V 3 ) —>■ { v r ,vp ,ip) where we have defined: tq =: v r , 
1> 2 =: vt cos {ip), vs =: vp sin {ip). It is easy to check that 
the Jacobian is given by: 


d {vi,V 2 , V 3 ) 
d {vi,v T ,ip) 


(Bl) 


The probability distribution f pk {v, v) is normalized, so the 
numerator in Eq. 16 can be expressed as: 


/ d 3 vf pk {v | v) = 1 - / d 3 vf pk {v | v) 

and we will now give a formula for the integral in the 
right-hand side of this equation. 

With the coordinate transformation given above we find, 
after some simple algebra: 


Md 


d 3 vf pk {v | v) = 


(a/47r 2 ) 


cO) 


a=i ^ 


1/2 


cO) 


dv r exp 


(v r - 


rq/^esc { r )— v r 


dvpvp exp 


dip exp 


1 1 
2 V^2 ^3 

1 1 


Ao A p 


{2iP) 


(B2) 


With the help of Eq. (9.6.16) from Abramowitz & Ste- 
gun (1972) the integral over the angular direction be¬ 
comes: 


dip exp 


1 1 


A •> A f 


K2 ip) 


= 2-irIn — 


1 1 

A2 A3 


(B3) 


where Jo denotes the modihed Bessel function of zero or¬ 
der. We now define the function: 


i-v 2 

F{r,\ 1 ,!3)= dze-P’Iofrz) 
Jo 

so that Eq. B2 becomes: 

/ d 3 vf pk {v,v) = —^-- 

J|t)|<t) eac 47T 11 j = 1 


(B4) 


c( r ) 


c( r ) 


exp 


(v r - . 


F (r] 2 {r),y,f]) (B5) 


d 2 (r) = v 2 esc {r)-v^, 


(B6) 


2 C 2 A 3 


P = 


1 f 1 


2 VA 


1 

A 3 


(B7) 


After having evaluated numerically the 
function F{r)2,y,l3) we can perform numerically the in¬ 
tegral (Eq. B5). This allows us to calculate the fraction of 
unbound particles. 
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